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ABSTRACT 

In this paper we study the effects of rotation on the frequencies of the quasi- 
normal modes of a proto-neutron star (PNS) born in a gravitational collapse 
during the first minute of life. Our analysis starts a few tenths of seconds after 
the PNS formation, when the stellar evolution can be described by a sequence 
of equilibrium configurations. We use the evolutionary models developed by 
Pons et al. (1999; 2001) that describe how a non rotating star cools down and 
contracts while neutrino diffusion and thermalization processes dominate the 
stellar dynamics. For assigned values of the evolution time, we set the star 
into slow rotation and integrate the equations of stellar perturbations in the 
Cowling approximation, both in the time domain and in the frequency domain, 
to find the quasi-normal mode frequencies. We study the secular instability of 
the g-modes, that are present in the oscillation spectrum due to the intense 
entropy and composition gradients that develop in the stellar interior, and 
we provide an estimate of the growth time of the unstable modes based on a 
post-Newtonian formula. 

Key words: gravitational waves — stars: oscillations — starsmeutron - 
stars: rotation — relativity - methods: numerical - 

1 INTRODUCTION 

It is well known that the frequencies of quasi normal modes (QNMs) of a star depend on 
its internal structure, and therefore on the particular evolutionary phase the star is going 
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through. In a recent paper (Ferrari, Miniutti & Pons 2003), to be referred to hereafter as 
FMP, it has been shown that these frequencies change during the first minute of life of a 
proto-neutron star (PNS) born in a gravitational collapse, and that the changes are mainly 
due to neutrino diffusion and thermalisation processes which smooth out the high entropy 
gradients that develop in the stellar interior. The models used in FMP were developed in 
Pons et al. (1999;2001) and describe the stellar evolution in terms of a sequence of equilibrium 
configurations; this 'quasi-stationary' approach has been shown to become appropriate a few 
tenths of seconds after the bounce following stellar core collapse which gives birth to the PNS. 
The study carried out in FMP shows that the frequency of all QNMs of a newly born PNS 
are much smaller than those of the cold NS which forms at the end of the evolution. Indeed, 
they initially cluster in a narrow region [y e [600, 1500] Hz) and begin to differentiate after 
less than a second. Unlike in zero temperature, chemically homogeneous stars, the frequency 
of the fundamental mode (f-mode) of a PNS does not scale as the square root of the average 
density, though the star is cooling and contracting; in addition, due to the strong thermal 
gradients that caracterize the initial life of the PNS, gravity modes (g-modes) are also 
present in the oscillation spectrum, and their frequencies are much higher than those of the 
core g-modes of cold NS's (Reisenegger & Goldreich 1992; Lai 1999). 

In FMP all stellar models were assumed to be non rotating; this choice was motivated by 
the need of isolating the effects of thermal and chemical evolution on the QNM spectrum. 
However, NSs are expected to be born with a significant amount of angular momentum, and 
the aim of this paper is to investigate how rotation modifies the picture described above. 
In particular, we shall explore the possibility that the g-modes become unstable due to 
Chandrasekhar-Friedman-Schutz (CFS) instability. This instability was first discovered by 
Chandrasekhar (1970) for the m = 2 bar mode of an incompressible Maclaurin spheroid, 
and was later shown to act in every rotating star by Friedman and Schutz (1978a,b). 

The CFS instability of the fundamental mode (f-mode) has been studied extensively in 
the literature in the framework of the newtonian theory of stellar perturbations (Bardeen 
et. al. 1977; Clement 1979), later generalized to the post-newtonian approximation (Cutler 
& Lindblom 1992), and more recently for fully relativistic, fastly rotating stars (Yoshida & 
Eriguchi 1997; Yoshida & Rezolla 2002). These studies show that the f-mode instability acts 
at very high stellar rotation rates, comparable to the break-up velocity limit of the star. It was 
also found that, unless the temperature is very low, viscous dissipation mechanisms tend to 
stabilise the f-mode instability. It should be stressed that, except that in Morsink, Stergioulas 
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& Blattning (1999) where a more realistic equation of state (EOS) has been considered, 
all these studies use a one parameter, polytropic EOS. The r-mode instability has also 
extensively been studied in recent years after Andersson (1998) pointed out that it is generic 
for every rotating star, and since the coupling of the r-mode with the current multipoles 
is strong, it was proposed that this instability plays an important role in nascent neutron 
stars. In principle, all quasinormal modes of a rotating star are unstable for appropriate 
values of the stellar angular velocity Q, and of the spherical harmonic index m, but the 
more interesting ones are those for which the instability sets in for a small value of Q. In this 
respect, it should be noted that since in the no rotation limit the g-modes have frequencies 
lower than the f-mode frequency, they may become unstable for relatively small values of 
the angular velocity. The g-mode instability has been studied only for zero temperature 
stars using the newtonian theory of stellar perturbations in the Cowling approximation. 
This study regarded a particular class of g-modes, for which the buoyancy is provided by 
the gradient of proton to neutron ratio in the interior of the star (Lai 1999). In this paper 
we study the onset of the CFS instability of the lowest g-modes of a newly born, hot 
proto-neutron star using, as mentioned above, evolutionary models that take into account 
the physical processes occurring in the early life of the star. We use the relativistic theory 
of stellar perturbations for a slowly rotating star in the Cowling approximation, which is 
known to reproduce with a good accuracy the g-mode frequencies because the gravitational 
perturbation induced by g-modes is much smaller than that associated to the f-mode. 

The paper is structured as follows. In section 2 we write the equations that describe a 
perturbed, slowly rotating, relativistic star both in the time and in the frequency domain in 
the Cowling approximation. In section 3 we present and discuss the numerical results both 
for the fundamental mode and for the lowest g-modes; we estimate the growth time of the 
unstable g-modes using post-newtonian formulae and draw the conclusions of our study. 

2 FORMULATION OF THE PROBLEM 
2.1 The background model 

We consider a relativistic star in uniform rotation with an angular velocity Q so slow that 
the distortion of its figure from spherical symmetry is of order Q 2 , and can be ignored. We 
expand all equations with respect the parameter e = Q/Qk, where VIk = djg, an d retain 
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only first order terms 0(e). On these assumptions, the metric can be written as (Hartle 
1976) 

ds > = _ e 2»(r) dt 2 + e 2A(r) rfr 2 + f .2 ^2 + sin 2 _ ^2 ^ 

The function u>(r) satisfies a second order linear equation 
4 / 4 \ 

ro ' ,r,r H — (A + l/) jr ( G7 jr H Wj = 0, (2) 

where we have defined 

w = Q — oo(r). (3) 

In the vacuum outside the star, A and v reduce to the Schwarzschild functions, and the 
solution of eq. (2) can be written as 

w = Q - 2Jr' 3 , (4) 

where J is the angular momentum of the star. The star is assumed to be composed by a 
perfect fluid, whose energy momentum tensor is 

Tp, = (p + p) u^u u + pg^, (5) 

with pressure p, energy density p and four- velocity components uF 1 = [e~ u , 0, 0, VLe~ u }. The 
metric functions u(r),X(r) are found by solving Einstein's equations for a spherically sym- 
metric, non rotating star, which couple the metric components to the fluid variables. As 
already mentioned in the introduction we shall use as a background the models of evolving 
proto-neutron stars developed in Pons et al. (1999;2001) and used in FMP; we shall choose 
different values of the evolution time starting from t = 0.5 s after the formation of the 
proto-neutron star, when the quasi-stationary description is appropriate to represent the 
stellar evolution. 

2.2 The perturbed equations in the Cowling approximation 

In this section we shall briefly outline the equations that describe the perturbations of a 
slowly rotating star up to first order in the rotation parameter e. We shall write these 
equations both in the time domain and in the frequency domain, because we have used both 
approaches to find the mode frequencies. In both cases we shall assume the oscillations to 
be adiabatic, so that the relation between the Eulerian perturbation of the pressure, 8p, and 
of the energy density, 5p, is given by 

5p=^5 P + p'e( I ^-l), (6) 

p + p v r / 
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where £ r is the radial component of the Lagrangian displacement and 1^ and T are 



p + p ( dp" 



r = 



p + p p' 



(7) 



P \ d PJ s,Y L 1 P P' ' 

a prime indicates differentiation with respect to r, and Yl is the lepton fraction. In the 
following equations we shall also use the speed of sound C s given by 

The complete set of the perturbed Einstein equations has been derived using the BCL 
gauge (Battiston, Cazzola & Lucaroni 1971) in Ruoff, Stavridis & Kokkotas (2002). The 
Cowling limit of these equations was studied in Ruoff, Stavridis & Kokkotas (2003), (RSK) 
for polytropic relativistic equations of state. In this work we will use the Cowling approx- 
imation, i.e. we shall neglect the contribution of the gravitational perturbations. In this 
approach, we need to consider only the fluid perturbations, i.e. the three components of the 
velocity perturbations Sui, the perturbation of the energy density 8p and the radial compo- 
nent of the displacement vector The perturbation of the pressure is related to 5p through 
the adiabatic condition (6). We expand the fluid perturbations as 



SUA; 



Lm 



(9) 



Lm 



f, lm f) V ~lm $ lm 

u 2 o e Yi m - u 3 . 

sin 



-e v ]T [u l ?d^Y lm + u l ™ sin 9d e Y lm 



l,m 



sp = E 



Lm 



(P + PY 



Y 



Imi 



V< 1- 



ri 



l,m 

where H = Sp/ (p + p) is the enthalpy. 



2.2. 1 Perturbed equations in the time domain 

The final set of perturbed equations in the time domain is 



(d t + imQ) H 



u 



, , 2;/ - A' + j) Ul - e 2A 4^ 2 + 2i///:,v 2A lv 11 



(d t + imQ) Ul = H' + 



P 



£i 
r 



1)H + Z 



Tip 

(d t + imft) u 2 = H + 2 ~^ ( imM 2 + Z^us) - ^-e -2 \Bui, 



- B (imu 2 + £f 1 u 3 ^j , 

2 

A 



(10) 
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2 

(d t + imQ) u 3 = 2^ (imu 3 - Cf 1 u 2 ) + ^e'^BC^m, 



where: 



A V * > K 

[1 

r 



A - 1(1 + 1), Q lm -\ {2l _ 1){2l + 1 y (11) 



b = J + 2w (V - ^ . 



In eqs. (10) we have omitted the indices Im in the perturbed variables and we have introduced 
the new quantities 

u 2 := u 2 + im— —e H, 



wr -21/ r ±i ti 
u 3 := u 3 — e £2 



2 

A" 

The operators Cf 1 , Cf 1 , and £3 1 are the same as in RSK and are defined by their action 
on a perturbation variable P lm 

C ±lplm = (|_l)Q k pW™_(I + 2 )Q i+lm pW™ (12) 

C ±lplm = _(I + l)g ira pi-l™ + I ft+lm pHlm 

C ±lplm = {l _ m+l)QlmP l-l m + l{l + 2)Qi+im pl + lm^ 

We must stress here that the above system is an infinite coupled system of differential 
equations ranging from I — m to I — 00 and that, in order to solve it numerically, we need 
to truncate the couplings to a finite value of Z max . 

In this work we want to study how the f-and g-modes, that have polar parity in the non 
rotating limit, are affected by rotation. When the star rotates, the equations that describe 
a polar perturbation with harmonic index / acquire rotational corrections with polar parity 
and index /, (for instance the terms imQH and 2imwe 2X ~ 2v H in the first of eqs. (10)) and 
with axial parity and index I ± 1 (for instance the term Cf l u 3 in the second of eqs. (10)). 

Since it has been shown that the polar rotational corrections are dominant for the modes 
under investigation (Kojima 1997), in our study we shall neglect the ±1 axial rotational 
corrections. Thus, the modes we are considering are those indicated in (Lockitch, Andersson 
& Friedman2001) as the "polar led hybrid modes". We choose / max = 2 and m = 2. We have 
checked that considering higher values of Z max would change the results by less than 1%. 

Eqs. (10) have been numerically integrated giving an initial gaussian pulse at the enthalpy 
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variable H . The frequencies of the modes are identified by looking at the peaks of the fast 
fourier transfom (FFT) of the outcoming signal. 



2.2.2 Perturbed equations in the frequency domain 

By replacing in Eqs. (10) all time derivatives by \a and letting H — > iH, w 3 — > fu 3 and 
£r _^ j£r^ we eas jiy obtain the real valued set of equations describing the eigenvalue problem 
in the frequency domain. We have two ODEs for H and u\ and three algebraic relations for 
u 2) u 3 and £ r . The relation for £ r , which follows from the last of eqs. (10) is particularly 
simple and can be used to eliminate that variable from the system. The result is 
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(a + mVL) U\ 



p 

Tip 
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(13) 



r 
A 



where we have defined 
2mm 



S := a + mil — 



A 



(14) 



An inspection of eqs. (13) shows that they become singular when E := 0. For any assigned 
value of /, m, Q and a, this may happen inside the star in a certain domain of the radial 
coordinate r which would depend on the values of the function m. This occurrence would 
generate the so-called continuous spectrum. As explained in RSK, for fixed values of the 
compactness of the star M/R, the frequency region were the continuous spectrum extends 
pratically depends on the values of m at the center and at the surface of the star and on the 
number of maximum couplings Z max that is considered. By following the procedure explained 
in RSK (section 2.3), it is easy to show that in the non-axisymmetric case (m ^ 0) and for 
Im&x = 2 the continuous spectrum extends to the following region 



m c 



2[ — -SI) < a < 2 — -Q 



m x 



H' 



it-, 



Assuming I — m — 2, Eqs. (13) can be written in the following simplified form 



2i/ + x'-- + c; 2 u' 
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(15) 



U! (16) 



mB 



Mi + 



2\-2u 



[2mm - C s 7 2 (o- + mQ)) + 



Ae 



2A n 



r 2 S 



H 



© 2002 RAS, MNRAS 000, 1-11 



8 V. Ferrari, L.Gualtieri, J. A. Pons, A. Stavridis 

where we have used the third of eqs. (13) to eliminate u-i- To find the mode frequencies, we 
integrate these two ODEs by imposing that the variables have a regular behaviour near the 
center, i.e. 

H ~ r\ Ul ~ r l -\ (17) 

and select those frequencies for which the Lagrangian perturbation of the pressure vanishes 
at the surface, i.e. 

Ap = (a + mVt)H(R) - u'(R)e 4u ^u 1 (R) = 0. (18) 



3 RESULTS 

In order to find the frequency of the quasi normal modes, we have numerically integrated 
the perturbed equations both in the time and in the frequency domain for different values 
of the evolution time t ev , and for selected values of the rotation parameter e. We choose the 
rotation rate to vary within < e < 0.4 because from preliminary calculations we find that 
for the models under consideration the mass shedding limit does not exceed e = 0.4 — 0.5. 

We consider the evolutionary model labelled as model A in FMP, in which the equation of 
state of baryonic matter is a finite-temperature, field-theoretical model solved at the mean 
field level. Electrons and muons are included in the models as non interacting particles, 
being the contribution due to their interactions much smaller than that of the free Fermi 
gas, and neutrino transport is treated using the diffusion approximation. The evolution 
time interval we consider covers the first minute of life of the proto-neutron star, from 
t ev = 0.5 s to t ev = 40 s, when processes related to neutrino diffusion and thermalization 
become negligible. The gravitational mass of the star at t ev = 0.2 s is M = 1.56 M , and at 
t ev = 40 s becomes M = 1.46 M Q . The difference in gravitational mass between the initial 
and final configuration is radiated away by neutrinos during the PNS evolution. The radius 
of the initial configuration is R — 23.7 km and reduces to R — 12.8 km at t ev = 40 s. 

For t ev <20 s the stellar models are convectively unstable, and the code which integrates 
the perturbed equations in the time domain explodes after some time, which is too short 
to accurately calculate the low frequencies of the g-modes. Conversely, the code which 
integrates the equations in the frequency domain is well behaved even when convective 
instability is present, and therefore for t ev < 20 s we use the frequency domain approach. 
After that time both methods can be applied and the results agree better than 5%. 
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It is worth stressing here once more that the perturbed equations in the frequency do- 
main present a singular structure which makes impossible their numerical integration in 
the continuous spectrum region. However, for the stellar models we use and for the mode 
frequencies we are interested in, we find that the continuous spectrum lays in the negative 
frequency range. 

The main results of this work are summarized in figures 1 and 2, where we plot the 
frequencies of the f-, gi- and g2- modes as a function of the rotation parameter e = Q/Qk, 
for different values of the evolution time in the more interesting phases of the cooling process. 

It should be reminded that the onset of the CFS instability is signaled by the vanishing 
of the mode frequency for some value of the angular velocity (neutral point). From figure 1 
and 2 we see that while the f-mode does not become unstable during the first minute of the 
PNS life, both the gi- and the g 2 - modes do become unstable. The gi- frequency remains 
positive during the first second, but at later times vanishes for very low values of e. For 
instance, at t ev = 3 s it crosses the zero axis for Q = 0.17 Qk, even though its value for 
the corresponding nonrotating star is still quite high, v gi = 486 Hz. The behaviour of the 
g 2 - mode is similar, but being the frequency lower the instability sets in at lower rotation 
rates. At later times, the g-modes frequencies decrease, reach a minimum for t ev = 12 s 
and then slightly increase. This behaviour can be attributed to the fact that during the first 
10-12 seconds the dynamical evolution of the star is dominated by strong entropy gradients 
that progressively smooth out. After about 12 s the entropy has become nearly constant 
throughout the star and g-modes due to composition gradients take over. 

It should be mentioned that in FMP we studied also a second model of evolving proto- 
neutron star, labelled as model B. The main difference between the two models is that model 
A has an EOS softer than model B, and that at some point of the evolution a quark core 
forms in the interior of model B. We have integrated the perturbed equations (10) and (16) 
also for model B, finding results entirely similar to those described above for model A; this 
indicates that that the quark core that develops at some point of the evolution in model B 
does not affect the overall properties of the modes in a relevant way. 

3.1 The growth time of the unstable modes 

A mode instability is relevant if its growth time is sufficiently small with respect to the 
timescales typical of the stellar dynamics, i.e. if the instability has sufficient time to grow 
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before other processes damp it out or the structure of the evolving star changes. In this 
section we will give an "order of magnitude" estimate of the growth time of the g-modes for 
the hot proto-neutron stars under investigation. Following Lai (1999), we shall evaluate the 
mode energy in the rotating frame using the expression given in Friedman & Schutz (1978a) 

E = \J p (a + mnf f • | + - 5p* + (V • () f • (Vp - C 2 s Vp) d 3 x, (19) 

and compute the growth time associated to the dissipative process we are considering, i.e. 
the gravitational radiation reaction, using the expression 

1 1 dE 

r gr 2E dt 

It is useful to remind that the relation between the Lagrangian displacement and the four- 
velocity of a perturbed fluid element is 

6u k = i(a + mn)e- v £ k , k = 1,3. (21) 

Since we are working in the Cowling approximation, the term 5$ in eq. (19) will be neglected. 
The energy loss due to gravitational waves can be calculated from the multipole radiation 
formula of Lindblom et. al. (1998) given by 

(it) =-°(° + mn )T, N l° 2l (\ SD l™\ 2 + \tJlm\ 2 ), (22) 
\ / gr l>2 

where the coupling constant N t is given by 
AnG (Z + !)(/ + 2) 

1 c 2m/(/_ 1)[(2/ + 1)!!] 2 ' 1 ' 

and SD [m and 5 J\ m are the mass and current multipoles of the perturbed fluid. The current 
multipoles 5Ji m are associated with the axial spherical harmonics; since we are interested 
in the g-modes, for which the effect of the coupling between polar and axial perturbations 
is negligible (indeed we did not include it in the perturbed equations), we shall neglect the 
8 Ji m contribution to the gravitational luminosity. The mass multipoles can be evaluated 
from the following integral expression 

SD lm = J r l 5pY lm d 3 x. (24) 

The growth times of the unstable gi- modes shown in figure 2 are summarized in Table 1. 
From these results we see that the growth time appears to be orders of magnitude larger than 
the timescale on which the star evolves, which is of the order of tens of seconds. Although 
we are aware that the estimate based on the newtonian expressions (19) and (24) is a quite 
crude one, the growth time is so much larger than the evolutionary timescale that it is 
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reasonable to conclude that the CFS instability of the lowest g-mode is unlikely to play any 
relevant role in the early evolution of proto-neutron stars. Similar conclusions can be drawn 
for the fundamental mode and for higher order g-modes. 
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Table 1. Growth times for unstable g\ mode of Model A for 
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Figure 1. The frequency of the fundamental mode of the evolving proto- neutron star is plotted as a function of the rotational 
parameter e = Q/Qk, for assigned values of the time elapsed from the gravitational collapse. We see that as the time increases, 
the frequency increases and tends to that of the cold neutron star which forms at the end of the evolutionary process. Reminding 
that the onset of the CFS instability occurs when the mode frequency becomes zero, we see that the f-mode would become 
unstable only for extremely high values of the rotational parameter, as it is for cold stars. 



© 2002 RAS, MNRAS 000, 1-11 



Rotational effects on the oscillation frequencies of newly born proto-neutron stars 13 




-500 



■1000 









*^3s 








\10s 






40 s 






1 1 1 




12s 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

£ 



N 
I 




■1500 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

£ 



Figure 2. The frequency of the gi- and g2- modes of the proto-neutron star are plotted, as in figure 1, as functions of the 
rotational parameter for assigned values of the time elapsed from the gravitational collapse. Unlike the f— mode, as the time 
increases the frequency of the g-modes decreases, reaches a minimum at about t ev = 12 and then slightly increases (see text). 
We see that for both modes the CFS instability sets in at values of the rotational parameter much lower that that needed for 
the f— mode. 
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